Analysis of self— averaging properties in the transport 
of particles through random media 

J.M. Lopez, M.A. Rodriguez and L. Pesquera 

Departamento de Fisica Moderna, Universidad de Cantabria, and Instituto de Estudios Avanzados en Fisica Moderna y 
Biologia Molecular, CSIC-UC, Avenida Los Castros, E-39005 Santander, Spain. 



> 

o 
in 

(N 



> 

o 






o 



% 



We investigate self-averaging properties in the transport of 
particles through random media. We show rigorously that in 
the subdiffusive anomalous regime transport coefficients are 
not self-averaging quantities. These quantities are exactly 
calculated in the case of directed random walks. In the case 
of general symmetric random walks a perturbative analysis 
around the Effective Medium Approximation (EMA) is per- 
formed. 
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The analysis of transport of particles in random me- 
dia has interest in physical systems since the transport 
mechanism is on the basis of many physical phenomena 
|^|-Q . The effect of disorder in the behavior of such sys- 
tems can be normal, if only a quantitative change of the 
transport parameters occur, or anomalous when a quali- 
tative change is induced by the disorder. This anomalous 
behavior is of great interest in the physics of disordered 
media and has been observed in almost all kind of phys- 
ical phenomena, from electrical conductivity to thermal 
properties P-H]. There is no general theory for anoma- 
lous transport because in general the involved phenom- 
ena are not unique. An important class of anomalous 
behaviors are those due to restrictions in the motion of 
particles imposed by the disorder. If this restriction is 
strong the motion of the particle is subdiffusive and ,as 
a consequence, anomalies in the observed phenomena oc- 
cur. 

An important characteristic associated to the anoma- 
lous behavior is the sample to sample dependence of the 
measured quantities. Strong sample to sample fluctu- 
ations are observed in most cases of anomalous behav- 
ior Q. In a normal situation the transport coefhcients 
are usually sample independent, while for strong disor- 
der transport coefficients are supposed not to be self- 
averaging quantities. The existence of sample to sam- 
ple fluctuations is problematic either from experimental 
or theoretical points of view. On the one hand, experi- 
ments are usually performed over few samples. On the 
other hand systematic analysis of disorder effects have 
been usually based on the probability density of diffus- 
ing particles averaged over random media configurations. 
The study of systems without self-averaging properties 
implies the calculation of averaged products of proba- 
bilities or some equivalent function, which is a rather 



difficult task. Only a few works have been devoted to 
the analysis of self-averaging properties in special cases. 
In ||] a case of weak disorder was investigated, while Q 
and [0 deal with cases of directed random walks. In 
recent works sample to sample fluctuations of first pas- 
sage times in asymmetric random walks S have been 
investigated. Also, and related with sample to sample 
fluctuations, there are some works dealing with sample 
averaging of powers of the probability ||^ . 

In this Letter we perform a systematic analysis of the 
self-averaging properties of systems described by random 
walks (RW) on regular lattices with a random distribu- 
tion of transition rates. We introduce a method based on 
the renormalization of some coefficients of the evolution 
equations to obtain averaged products of probabilities. 
This method generalizes the one previously used in llG] 
for the calculation of single averaged probabilities. We 
apply the method, which we believe to be of rather wide 
applicability, to the general symmetric RW and to the di- 
rected RW in one dimension. It is shown rigorously that 
for strong quenched disorder the transport coefficients 
arc not self-averaging. 

We consider a general transport problem in which a 
particle moves in a lattice with random transition rates. 
The position of the particle at a time t is denoted by 
r{t) and the magnitudes of interest are mean functions 
of the position F(r). Defining P{r, t) as the probability of 
finding the particle in r at time t, the observed quantities 
are given by: 



F(i)-^F(r)P(r,i). 



(1) 



In principle these quantities are dependent of the particu- 
lar configuration of the medium. A complete description 
of these quantities can be achieved by using averaged 
moments of the form: 

(F"(t))- E i^(ri)...F(r„)(P(ri,t)...P(r„,t)), (2) 



where (...) indicate average over all possible configura- 
tions of the medium. The self-averaging character of F{t) 
can be derived from its variance. A zero dispersion is 
equivalent to self-averaging. This dispersion can be ob- 
tained from the calculation of the averaged products of 
probabilities. In the following we focus on this problem 
in a transport model with a probability governed by a 
master equation with random coefficients Wr as: 



dtP{R, t) = L^P + LLWrLnP, 



(3) 



where Lq, Ll and Ln are linear operators. In a stan- 
dard transport model these operators are linear combi- 
nations of shift operators, L{r) = '^iaiEi{r), such that 
Ei{r)P{r^t) = P{r + i,t) and Wr are random transition 
rates. 

The starting point of our method is the introduction 
of an effective medium with memory |IQ] . In this way (H) 
can be written in terms of Laplace Transforms as 

sP{r, s) - Po{r) = [Lo + Ll<P{s, r)LR]P{r, s) + L^K - </* 

(4) 

where Pq (r) is the probability at i = and (l){r, s) is the 
transition probability of the effective medium that will 
be determined below. Taking the integral form of (Q) 
and iterating, a development in powers of the random 
transition rate 9sir) = Wr — (f>{s, r) is obtained ||l^. Now 
we renormalize 9sir) by performing a summation of all 
terms in which contiguous indexes take the same value 
(|lO| , obtaining 

oo 

P{r, s) = Gs{r, r/)Po{r/) + ^ $,(ri) ^s{rn)G^{r, n)J, 

n=l 

(5) 

where summation over repeated indexes is understood 
and the sum is restricted to terms with different contigu- 
ous indexes. This renormalization corresponds to the one 
loop resummation in diagrammatic representations, also 
known as single site approximation in condensed matter. 
The renormalized random transition is given by 



difference between the exact probability and that ob- 
tained with the effective medium: SP{r, s) = P{r, s) — 
Gs{r,r/)Po{r/). In this way the averaged products of 
SP{r, t) are obtained from (g) as series in moments of 
^^(r). Since self-averaging is equivalent to a null disper- 
sion, to analyze the self-averaging character of the trans- 
port coefficients only {P{r,s)) and {P{r,s)P{rf,s)) must 
be considered. The method outlined above can be used 
to obtain these magnitudes in a large variety of prob- 
lems. Here we consider the directed RW and the general 
, symmetric RW in one-dimensional media with quenched 

a. Directed random walk (DRW) in ID. In the DRW 
only steps in one direction are allowed. Despite its sim- 
plicity several phases or anomalous behaviors appear de- 
pending on the intensity of disorder |6| . In one dimension 
the master equation modeling the DRW can be written 
as: 



dtPin,t) = -(1 - E^i{n))wnPin,t). 



(9) 



The anomalous phases can be classified according to 
the intensity of disorder, which is related to the existence 
of inverse moments ofjthe random term w„ . If we restrict 
(^t(ir]?W;il5/fe&-tb't%l'^li^trnil'oi^il^avior of the velocity, 
only the existence of the first inverse moment is relevant. 
Taking a probability distribution p(w„) = (1 — a)w~°' 
the weak disordered phase correspond to the existence 
of the first inverse moment, (a < 0), and the strong 
disordered phase to {w~^) ^ cx), (1 > a > 0). Other 
cases concerning transients can be found in H . 

The application of the method to this case is straight- 
forward. The propagator Gs{n,m) is zero when n < m 
and for n > m we have 



^s{r) 



Os{r) 



1-Js{r,r)es{r) 
and the functions G^'^ and Jg are defined by 

Gf{r,rf) = LR{r)Gs{r,rf),G^{r,r/) = Ll{r/)Gs{r,r/) 



(6) 



Gs{n,m) 



(s + 0(3))"-™+! ■ 



(10) 



Mr,r/)^Ln{r)Ll{rf)Gs{r,rf), 



(7) 



(8) 



being G's(r, r') the propagator of the deterministic part 
of (^, Lq + LL(j)Lpi. Finally, the transition probability 
(j){r, s) is defined by the Effective Medium Approximation 
(EMA) condition [|o|: ($s(r)) = 0. When the model is 
translationally invariant the propagator is only depen- 
dent on the difference of site positions and the effective 
medium is homogeneous, that is, (/> is not dependent on 
the position. 

The averaged products of probabilities can be directly 
calculated from (ph. These products are more conve- 
niently expressed in terms of 6P{r,s), defined as the 



Using the EMA condition we obtain the transition prob- 
ability of the effective medium (j>{s) = R~^{s) — s where 
the function R{s) = {{s + w)~^) has been calculated 
in Ref. H. Since in the DRW only steps in one di- 
rection are possible, only terms in (0) with ordered in- 
dexes ri > r2 > ... > r-n are different from zero. Then 
{SP{r, s)) = and {P{r, s)) is exactly given by the EMA. 
It is also possible to obtain exact expressions for the av- 
eraged products of the moment generating function de- 
fined as F{x, s) = X^i^o ^^P{h s)- The factorial moments 
/„(s) can be obtained by taking the derivative of F{x, s) 
at a; = 1. All these quantities are sample dependent. The 
averaged products of factorial moments can be calculated 
by means of averaged products of generating functions as 



(/«(s).../m(s)) = 



a" 



'(F(x,s)...F(y,s)) 



dx''K..dy" 



(11) 



In general any self-averaging property can be analyzed 
with the knowledge of the averaged products of generat- 
ing functions. Let us consider as an example the analysis 
of the behavior of the factorial moments and assume the 
asymptotic form fn{s) ~ Uns"'" , where a„ is in principle a 
sample dependent quantity. The averaged value of a„ and 
its dispersion can be obtained from {F(x, s)) and from 
D{x,y,s) = {{Fix,s) - {Fix,smFiy,s) - {Fiy.s)))). 
From (ph we obtain the exact expression of the averaged 
products of generating functions: 



D(x,y,s) 



{F{x,s)) 



($? 



1 



0(s)(l-x) 



{l-x){l-y) 



sample averaged magnitudes are the same for RT and RB 
in one dimension [nOl . The same results are also obtained 
from (5) for the self-averaging properties. 

The application of the method to RB and RT models 
is also straightforward, but it is not possible to obtain 
exact expressions like in the DRW case. The propagator 
Gs{n, m) and the functions Js{n, m) are given in llG] for 
all kinds of disorder. The transition probability of the ef- 
fective medium </)(s) has been also calculated in ||l^] from 
the EMA condition. In this reference we obtained the 
exact asymptotic behavior of the averaged mean square 
(12) displacement (a;^(s)) in the frequency domain, which is 
directly related with the diffusivity. The results derived 
from the EMA for each type of disorder (A, B and C) are 
give$ by : 



(S + 0(S))2 [.S+{l-x)(b{.s)][s+il 



(13) 



(1 



where ^(s) = [0(s)2(s + 0(s)) s^^*^^] (s + 0(s)) 
and the renormalized random transition is <I>s = (s -|- 
(f>{s))[(w — (/)(s))/(w -|- s)] . From these expressions it 
is immediate the calculation of moments and their sam- 
ple to sample dispersions. In the weak disordered phase, 
after calculation of R{s) and {^'^), one obtains a ballis- 
tic behavior, {r(t)) ^ vt, with a velocity v = {w~^)~^, 
that is a self-averaging quantity. In the strong disor- 
dered phase one obtains, in agreement with [6], [7], a 
subballistic behavior, {r{t)) '^ 6^^^^") with a coefficient 
with mean value (6) — sin(7r(l — a))/(7r(l — a)r{2 — a)), 
that is not self-averaging. The relative variance of b is 
a^{b) = {{b - {b)f)/{b)^ = a/{2 - a). The dispersion 
increases {ar{b) -^ 1) for stronger disorder {a — > 1). 

b. Symmetric random, walk (SRW) in ID. There are 
two models of symmetric RW, the random trap (RT) and 
random barrier (RB) models. The master equations cor- 
responding to both models are written in our formulation 
as: 



A{s)xy) ' 
)ema 



x^{s)i 



^Pi 



1 



2/3? 



(/3is)U0(s)) (16) 



x'^{s) 



EMA 



s^ I ln(s 



4 In I ln(,s) 



In; 



+ 0(1 In. 1-1)) 
(17) 



x'^{s)ema = 2 



sin(7ra) -12/(2-0!) 3e,_4 sc-e . 



(1 -a)7r2" 



(18) 



The exact results can be expressed as corrections to 
the results given by the EMA as: 



(a;2(s)) r^x'^{s)EMA{l + a{s)) 
where the first corrections are 



O-Ais) 



1 



12/3f 



(/32 - PiYs 



(19) 



(20) 



dtP[n,t) = (1 - E^i[n))w^{E+i{n) - l)P(n,t) (14) 

for the random barrier and 

dtP{n,t) = {E^i{n)+E+i{n)-2)wr,P{n,t) (15) 

for the random trap. The anomalous behavior induced 
by the disorder is, in both cases, well known |1|. The 
different phases can be classified following the definitions 
of [1]. We recall that for model A (weak disorder) the 
inverse moments of Wm Pm = {{wn)~^') {M = 1,2,...) 
arc finite, while models B (marginal case, a = 0) and 
C (strong disorder, < a < 1) are based on a proba- 
bility distribution p{wn) = (1 — a)w~°' (w„ G (0,1)), 
such that inverse moments diverge. In all cases the long 
time behavior of the sample averaged diffusion coefficient 
has been exactly calculated. However sample to sample 
fluctuations have not been investigated until now. The 



as(s) 



(tt^ + 16 In 2 - 20) 
llnsP 



acis) = (41n2-5 + 7rV4)a2 + 0(aMna). 



(21) 



(22) 



In the weak disordered case the behavior is normal and 
the EMA reproduces exactly the first and second terms 
of (^2). In the marginal case B the EMA is exact up 
to terms of order smaller than | Ins |~^ pM. In the 
strong disordered case the behavior is subdiffusive and 
the EMA does not reproduce exactly the coefficient of 
the leading term pn]. Expressions (19)-(22) have been 
diagrammatically calculated in [Q by using cumulants 
and projection operators. We can obtain the same re- 
sult from (^ in a much more simple way in terms of 
moments and single functions. This simplicity allows us 



to calculate more involved quantities and to analize self- 
averaging properties. For instance to analize the sample 
to sample fluctuations of the generalized diffusion coef- 

2 

ficient we have calculated {x^ ) which depends on the 
averaged product of probabilities. As in the above case 
the exact result can be expressed as corrections to the 
EMA result as 



where the correction terms are 



la 



-S2 



1 



lb 



\-D.i 



7c = a/2 + 0(a2). 



(23) 



(24) 



(25) 



(26) 



These terms have been calculated from a diagrammatic 
representation of the averaged products of (H) . A detailed 
description of this method will be presented elsewhere. 
Finally, from (nS) and (ESI) we can extract several con- 
clusions. For weak and strong disorder the dispersion of 
the particle can be taken, in the long time limit, as x'^ ^ 
ais"^ -|-a2s"^ , where the coefhcients are in principle sam- 
ple dependent quantities. In the weak disordered case the 
behavior is diffusive and we have x'^ ~ ais~^ -\- a-is '^ ^ 
where the first coefficient is self-averaging but the second 
is sample dependent with a zero mean value and a vari- 
ance 0^(0,2) = Pi {(32 — 0i)- In the strong disordered 
case the behavior is subdiffusive, x"^ ~ cis~('*~'^"^/(^~") 
and ci is not self-averaging. The mean value of the co- 
efficient ci can be easily calculated for small a from (|l9| ) 
and (22) ((ci) = 4a -I- 0(a^)). As we will show else- 
where its dispersion can be also obtained in the same way 
cr^(ci) — 2a^ ~\-0{a^ Ina). Finally, the marginal case B is 
similar to the weak case with logarithmic corrections and 
we have x^ ~ 6is~^ | Ins |~^ +123^"^ \ Ins |~^/^. The 
first coefficient is self-averaging but 62 is sample depen- 
dent with a zero mean value and a variance (7^(62) = 16. 

In summary, we have presented here a general method 
to study self-averaging properties in the transport on 
random media. Our analysis of both DRW and SRW 
shows rigorously that when the behavior of a magnitude 
is normal its long time behavior is sample independent. 
By the contrary in anomalous diffusion phases the self- 
averaging property is not satisfied. The method intro- 
duced in this letter can be easily applied to other one di- 
mensional problems and it can be also extended to more 
dimensions. Some of these applications will be presented 
elsewhere. 
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